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A quench in an underdamped one dimensional cj> 
model is studied by analytical methods. The density 
of kinks just after the transition is proportional to 
the square root of the rate of the quench for slow 
quenches. If the quench is shorter that the relaxation 
time, then the density scales like the third root of the 
rate. 

Topological defects play a prominent role in many con- 
densed matter systems, see e.g. jl| for a review. It was 
also suggested that they were an important ingredient of 
the early universe 0. Topological defects can be gener- 
ated in large numbers during a second order phase tran- 
sition. The dynamics of such a transition to a symmetry 
broken phase has been an object of much recent attention 
because of its importance in the cosmological context || 
and in condensed matter physics B . An early estimate of 
defect density after a quench was given by Kibble [9. It 
his theory, the speed of light is a dominant factor which 
determines the size of correlated domains. Another sce- 
nario was put forward by Zurek || , who emphasized the 
importance of the noncquilibrium dynamics of the order 
parameter. Recent experiments point to the latter the- 
ory. 

The key observation we make in this communication 
is that the relaxation time in the underdamped limit is 
the same for all Fourier modes and it is an inverse of the 
viscosity coefficient 7. If a quench is slow enough so that 
the system can equilibrate during the initial stage of the 
quench, then this relaxation time determines the frozen 
correlation length and thus the density of kinks just af- 
ter a quench. If a quench is faster than the relaxation 
time, then there is not enough time to equilibrate and 
the system developes a characteristic length scale out of 
equlibrium. 

To be more specific, let us consider a one dimensional 
underdamped </> 4 model 



4>(t,x) +7^0, x) = 
(f>"{t,x) - 2a(t)<t>(t,x) 



2((> 3 (t,x) +rj(t, x) 



(1) 



where ' = dt, ' = d x . <f>(t,x) is a real order parameter. 

The system is assumed to be strongly underdamped for 

2 

the relevant part of the quench, 2a(t) » \. r](t,x) is a 
Gaussian noise of temperature T with correlations 



< rj(t,x) >= , 

< rj(ti,xi)r](t2,X2) >= 2jTS(ti 



t 2 )6{x 1 -x 2 ) • (2) 



The coefficient a(t) is time dependent. We consider a 
linear quench 



a(t) 




if t < 

if < t < r^±i 

if T±£ 



(3) 



< t 



Before the quench, for t < 0, the system is in a symmetric 
phase (A > 0), during the quench, at t — t, it undergoes 
a transition from the symmetric phase a{t < r) > to a 
broken symmetry phase a(t > t) < 0. Finally it settles 
down at a(t) = —1. 

As long as the system is still in the symmetric phase 
(t < t) and the temperature is moderate, the field </> can 
be regarded as a small fluctuation around its symmetric 
ground state < (j)(t,x) >= 0. It is justified to neglect 
the cubic term on the RHS of Eq.(§). The field can be 
written as a Fourier transform, 

/+00 
dk^j)(t,k)e lkx . (4) 
-00 

The Fourier transform of linearized Eq.(^j) and of the 
noises (||) is 

4>(t, k) + j4>(t, k) = ~[2a(t) + k 2 ]4>{t, k) + i%t, k) , 

< f]{t,k) >= , 

< fj*{h, fci)jj(i 2l k 2 ) >= ^-5{h ~ t 2 )5{k x - k 2 ) . (5) 

To predict the density of kinks in the symmetric phase 
we must know the two point correlation function < 
<p*(T,ki)<p(T,k 2 ) > at the critical point t — t. An ex- 
act formal expression can be constructed for any 7 as a 
rather involved combination of Bessel functions but it is 
too complicated to be very illuminating. We concentrate 
here on the limit of very small 7, where some intuitive 
understanding can be achieved. 

The limit 7 — > is a delicate matter as 7 appears 
both in the field equations and in the noise correlations. 
Any ad hoc setting of 7 = here or there may be dan- 
gerous. Our approach to this problem is as follows. In 
the strongly underdamped limit the system (|^) can be 
mapped to the equivalent one 
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4>(t, k) 



n(t, k) 



1 



< n(t, fc) >= , 

< n{t 1 ,ki)n(t 2 ,k 2 ) >= 



~fT5(h - t 2 )S(k 1 ~ k 2 ) 



fc 2 ] 



(6) 



2n[2a(t 1 

If a(t) = a = const, then both systems have correlations 

r<5(fci - hi) 



<0*(t,fclMt, fc 2 ) >= 



27r[2a + fc 2 ] 



(7) 



Their response to the change of parameters is the same 
too. In particular they have the same relaxation time. 
For example, if a(t) = a = const and the temperature 
jumps from T = for t < to T > for t > 0, then the 
correlation functions for t > are 

< 4>*(t, kl )4>(t, k 2 ) >= ggg 7 §j [1 - e-^l 



2?r[2a + fc 2 ] 



< k x )4>{t,k 2 ) > 
T8(ki - ha) 



1 



,-7< 



0(- 



7 



r)] 



(8) 



27r[2a + fc 2 ] L ~ V2a + fc 2 

for the models (^|) and (JsJ) respectively. 0( , 2 J ) de- 
notes oscillating terms which are negligible for small 7. 
In this limit the response of the two systems is the same. 
Another example is the response to a sudden change in 
the parameter a from a = a_ for t < to a = a + for 
i > 



<^*(t,fciMi,fe 2 ) >= 



T5(fci - fc 2 ) 1 - e 



- 7 t 



- 7 i 



2tt 



L 2a 1 



A: 2 2a_ 



< 4>*(t,h)4>(t,h) >= T6{h k2) x 



- 7 t 



-7* 



L 2a+ + k 2 2a_ + fc 2 



2tt 
0( 



7e 



7e' 



-7T 



^J2a + + k 2 y/2a- + h 



where the 0(. . .) terms are once again oscillations negligi- 
ble for 7^0. The systems respond in the same way to an 
impulsive perturbation of their parameters. The impul- 
sive perturbation contains all frequencies so the models 
are equivalent in the underdamped limit. 

From now on we work in the framework of the model 
(p|). The formal solution of Eqs.dq) is 



0(r, fc) 



dt\ e 



- 7 (-r-ti)/2 T 



ih,h) 



(10) 



This solution, the correlations (||) and the explicit form 
of alt), see Eq.(^), give the correlation function at t = r 



< 4>* (r, kt)4>( T , fc 2 ) >= ^S(hx - fc 2 ) G(r, fcO 



G(r, fc) 



-7T 



2A + k 2 



dx 



^x + k 2 

7T 



(11) 



G(t, fc) has a very interesting asymptote for r >> 1/7 

1 



G(r,k) 



2A 
7T 



fc 2 



(12) 



This asymptote can be understood in heuristic way. In 
the extreme underdamped limit (^) all the Fourier modes 
have the same relaxation time 1 /j whatever is their mo- 
mentum fc. After the quench begins at t — the system 
tries to adapt to the changes of the parameter a(t). It can 
not react on a time scale shorter than 1/7. The correla- 
tions of the system are effectively " frozen" at the instant 
t, when the time still left to the transition at t = r is 
equal to the relaxation time, (r — t) — I/7. At this in- 
stant a(i) = — and the correlation function is frozen as 

V ' 7T 

(|l2|). The frozen correlation length is so tne density 
of kinks just after a quench scales like 



(13) 



provided that the quench time is longer than the relax- 
ation time 



fc 2J 



:)] , ( 




T » 



1 

1 



(14) 



and it is short enough so that the system is still under- 
damped at the freeze-in time t, 2a(t) » 7 2 /4 or 



t << 



8A 

^y3 



(15) 



The range of quench times r where the t > scaling 
applies blows up as 7 tends to zero. 

If the quench time is shorter than the relaxation time 
t << I/7 then the system is out of equlibrium from 
the very beginning of the quench. As the system does 
jaot equilibrate anyway, we can effectively set 7 — in 
Eqs.(y,||). In the absence of damping exact solutions of 
equations (^]) for the linear quench @ are combinations 
of the modes 



4> k {o < t < t) 



-(1 



fc 2 ^ 

2A' 



-t J 




21 (r(l + |[)-t)«] . (16) 



At the critical point t = t the modes take the form 



rfc 2 r T ^ 3 i 
~2A ±i [ 3A J 



(17) 



The correlation function a,t t = r can be constructed as 
a combination of cf>k(0), 4> k (0) and of <^fc(r). The latter 
contain information about a characteristic length scale 
at the critical point. At A- w 1 the behaviour of the 
Bessel function (|l7|) changes qualitatively. An inverse of 
this length scale determines the density of kinks just after 
a quench 



(-)' 

T 



(18) 



provided that 



2 



1 , X 

T « - . (19) 

7 

This result is consistent with the discussion of Laguna 
and Zurek and with results of Karra and Rivers || . 
The opposite r >> 1/7 regime has not been discussed 
before. 



Conclusion. 

For fast underdamped quenches the frozen correlation 
length grows like r 1 / 3 . For slow underdamped quenches 
it scales like r 1 / 2 . 
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